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ABSTRACT. Given an ensemble of N x N random matrices, a natural question to ask is whether 
or not the empirical spectral measures of typical matrices converge to a limiting spectral measure 
as N — > oo. While this has been proved for many thin patterned ensembles sitting inside all real 
symmetric matrices, frequently there is no nice closed form expression for the limiting measure. 
Further, current theorems provide few pictures of transitions between ensembles. We consider the 
ensemble of symmetric m-block circulant matrices with entries i.i.d.r.v. These matrices have toroidal 
diagonals periodic of period m. We view m as a "dial" we can "turn" from the thin ensemble of 
symmetric circulant matrices, whose limiting eigenvalue density is a Gaussian, to all real symmetric 
matrices, whose limiting eigenvalue density is a semi-circle. The limiting eigenvalue densities f m 
show a visually stunning convergence to the semi-circle as m —> oo, which we prove. 

In contrast to most studies of patterned matrix ensembles, our paper gives explicit closed form 
expressions for the densities. We prove that f m is the product of a Gaussian and a certain even 
polynomial of degree 2m — 2; the formula is the same as that for the m x m Gaussian Unitary En- 
semble (GUE). The proof is by derivation of the moments from the eigenvalue trace formula. The 
new feature, which allows us to obtain closed form expressions, is converting the central combina- 
torial problem in the moment calculation into an equivalent counting problem in algebraic topology. 
We end with a generalization of the m-block circulant pattern, dropping the assumption that the m 
random variables be distinct. We prove that the limiting spectral distribution exists and is determined 
by the pattern of the independent elements within an ?7i-period, depending on not only the frequency 
at which each element appears, but also the way the elements are arranged. 
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1. Introduction 

1.1. History and Ensembles. Random matrix theory is the study of properties of matrices chosen 
according to some notion of randomness, which can range from taking the structurally independent 
entries as independent identically distributed random variables to looking at subgroups of the classi- 
cal compact groups under Haar measure. While the origins of the subject go back to Wishart's HWisll 
investigations in statistics in the 1920s, it was Wigner's work [ Wigl[ Wig2[ Wig3[ Wig4[ Wig5[ in 



the 1950s and Dyson's [ Dyl[ Dy2[ a few years later that showed its incredible power and utility, as 



random matrix ensembles successfully modeled the difficult problem of the distribution of energy 
levels of heavy nuclei. The next milestone was twenty years later, when Montgomery and Dyson 
HMonll observed that the behavior of eigenvalues in certain random matrix ensembles correctly de- 
scribe the statistical behavior of the zeros of the Riemann zeta function. The subject continues to 
grow, with new applications ranging from chemistry to network theory HMNSM to transportation 
systems HBBDSl IKrSeL See HFM[ |Hay[ for a history of the development of the subject and the 
discovery of some of these connections. 

One of the most studied matrix ensembles is the ensemble of N x N real symmetric matrices. 
The N entries on the main diagonal and the ^N(N — 1) entries in the upper right are taken to 
be independent, identically distributed random variables from a fixed probability distribution with 
density p having mean 0, variance 1, and finite higher moments. The remaining entries are filled in 
so that the matrix is real symmetric. Thus 

[Pa 

Prob(A) = || P(aij), Prob(A : ay G [c%,Aj]) = / i>{.r,j)(I.r, r (1.1) 

l<i<j<N l<i<j<N J x ij= a ij 

We want to understand the eigenvalues of A as we average over the family. Let 5{x — x ) denote 
the shifted Delta functional (i.e., a unit point mass at x , satisfying J f(x)5(x — x )dx = f(x )). 
To each A we associate its empirical spacing measure: 

i=l ^ V / 
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Using the Central Limit Theorem, one readily sees that the correct scale to study the eigenvalues is 
on the order of v^vQ The most natural question to ask is: How many normalized eigenvalues of a 
'typical' matrix lie in a fixed interval as N — > oo? Wigner proved that the answer is the semi-circle. 
This means that as N — > oo the empirical spacing measures of almost all A converge to the density 
of the semi-ellipse (with our normalization), whose density is 



fwi S {x) = <( ^ V i i - (i.3) 

otherwise; 



to obtain the standard semi-circle law we need to normalize the eigenvalues by 2 v N and not y N. 

As the eigenvalues of any real symmetric matrix are real, we can ask whether or not a limiting 
distribution exists for the density of normalized eigenvalues for other ensembles. There are many 
interesting families to study. McKay HMcKI proved that the limiting spectral measure for adjacency 
matrices attached to d-regular graphs on N vertices exists, and as N — > oo, for almost all such 
graphs the associated measures converge to Kesten's measure 



,l^)V^d-l)-x 2 , \x\<2Vd^l 

jKcstcn,d{X) — < y ' . (.1.4) 

otherwise 

(note that the measures may be scaled such that as d — > oo they converge to the semi-circle distri- 
bution). 

This example and its behavior are typical for what we hope to find and prove. Specifically, we 
are looking for a thin subfamily that has different behavior but, as we fatten the ensemble to the full 
family of all real symmetric matrices, the limiting spectral measure converges to the semi-circle. 
Numerous researchers have studied a multitude of special, patterned matrices; we do not attempt to 
do this vast subject justice, but rather concentrate on a few ensembles closely related to our work. 

All of the ensembles we consider here are linked ensembles (see HBanBolO . A linked ensemble of 
N x N matrices is specified by a link function Ljy : {1,2,..., N} 2 — > S to some set S. To s G S, 
assign random variables x s which are independent, identically distributed from a fixed probability 
distribution with density p having mean 0, variance 1, and finite higher moments. Set the (i, j) th 
entry of the matrix a it j := xz^jjjH For some linked ensembles, including those we examine here, 
it is be more convenient to specify the ensemble not by the link function, but by the equivalence 
relation ~ it induces on{l,2,...,iV} 2 . A link function may be uncovered as the quotient map to the 
set of equivalence classes {1,2,..., N} 2 -» {1,2,..., N} 2 / ~. For example, the real symmetric 
ensemble is specified by the equivalence relation ~ (j, i). 

One interesting thin linked ensemble is that of real symmetric Toeplitz matrices, which are con- 
stant along its diagonals. The limiting measure is close to but not a Gaussian (see HBCGl IBDJl 
IHMIO ; however, in HMMSM the sub-ensemble where the first row is replaced with a palindrome is 
shown to have the Gaussian as its limiting measure. While the approach in HMMSI involves an 
analysis of an associated system of Diophantine equations, using Cauchy's interlacing property one 



^Y^iLi = Trace(A 2 ) = j<jy a ?j' as ^ e mean is zero and the variance is one for each a,ij, this sum is of the 
order TV 2 , implying the average square of an eigenvalue is N. 

For general linked ensembles, it may make more sense to weight the random variables by how often they occur in 
the matrix: tijj := ciq\L^{{Liq{i, j)})| _1 ^L JV (ij)- For the real symmetric ensemble, this corresponds to weighting 
the entries along the diagonal by 2. In that case, and for the ensembles we examine here, this modification changes only 
lower order terms in the calculations of the limiting spectral measure. 
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can show that this problem is equivalent to determining the limiting spectral measure of symmetric 
circulant matrices (also studied in HBMIO . 

While these and other ensembles related to circulant, Toeplitz, and patterned matrices are a very 
active area IBasBoll IBasBo2l IBanBol iBCGl IBH1 IBM iBDJl IHMllMMSL of particular interest to 
us are ensembles of patterned matrices with a variable parameter controlling the symmetry. We 
desire to deform a family of matrices, starting off with a highly structured family and ending with 
the essentially structureless case of real symmetric matrices. This is in contrast to some other work, 
such as Kargin UKarll (who studied banded Toeplitz matrices) and Jackson, Miller, and Pham HJMPH 
(who studied Toeplitz matrices whose first row had a fixed but arbitrarily number of palindromes). 
In these cases the ensembles are converging to the full Toeplitz ensemble (either as the band grows 
or the number of palindromes decreases). 

Our main ensemble is what we call the ensemble of m-block circulant matrices. A real sym- 
metric circulant matrix (also called a symmetric circulant matrix) is a real symmetric matrix that 
is constant along diagonals and has first row (x , X\, £2, • • • , ^2, Note that except for the main 
diagonal, a diagonal of length N — k in the upper right is paired with a diagonal of length k in 
the bottom left, and all entries along these two diagonals are equal. We study block Toeplitz and 
circulant matrices with m x m blocks. The diagonals of such matrices are periodic of period m. 

Definition 1.1 (m-Block Toeplitz and Circulant Matrices). Let m\N. An N x N real symmetric 
m-block Toeplitz matrix is a Toeplitz matrix of the form 



( B 


B l 


B 2 ■ 


■ B N / m -i 


\ 


B-i 


Bo 


Bi ■ 


BN/ m -2 




B-2 


B-! 


Bo ■ 


BN/ m -3 




\ Bl_N/ m 


B2-N/ m 




■ Bo 


J 



with each Bi anmxm real matrix. An m-block circulant matrix is one of the above form for which 

B—i -B n _j. 

We investigate real symmetric m-block Toeplitz and circulant matrices. In such matrices, a 
generic set of paired diagonals is composed of m independent entries, placed periodically; however, 
as the matrix is real symmetric, this condition occasionally forces additional entries on the paired 
diagonals of length N /2 to be equal. 

For example, an 8 x 8 symmetric 2-block Toeplitz matrix has the form 



/ C Ci 

ci do 


c 2 c 3 
di d 2 


C4 C5 
d 3 d 4 


c 6 c 7 \ 
d 5 do 


c 2 d 1 
c 3 d 2 


Co Ci 

Ci do 


c 2 c 3 
di d 2 


C4 C5 
d 3 d 4 


c 4 d 3 

C5 G?4 


c 2 di 
c 3 d 2 


Cq Ci 

ci d 


c 2 c 3 
di d 2 


c 6 d 5 
\ c 7 d 6 


c 4 d 3 
c 5 d 4 


c 2 di 
c 3 d 2 


Co Ci 

ci d J 



(1.5) 
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while a 6 x 6 and an 8 x 8 symmetric 2-block circulant matrix have the form 



/ c 


Cl 


c 2 


c 3 


c 2 


di 


\ 


Cl 


d 


d x 


d 2 


c 3 


d 2 


c 2 


di 


c 


Cl 


c 2 


c 3 




c 3 


d 2 


Cl 


do 


di 


d 2 




c 2 


c 3 


c 2 


di 


Co 


Cl 




V di 


d 2 


c 3 


d 2 


Cl 


d 


/ 



( C Ci 

Ci do 


c 2 c 3 
di d 2 


j 

c 4 « 3 
d 3 d 4 


c 2 di \ 
c 3 d 2 


c 2 di 
c 3 d 2 


Co Ci 

ci d 


c 2 c 3 
di d 2 


c 4 d 3 
d 3 d 4 


c 4 d 3 
d 3 d 4 


c 2 di 
c 3 d 2 


Co Ci 

ci d 


c 2 c 3 
di d 2 


c 2 c 3 
\ di d 2 


c 4 d 3 

G?4 


c 2 di 
c 3 d 2 


Co Ci 

ci d / 



(1.6) 



Note for the 6x6 matrix that being real symmetric forces the paired diagonals of length N / 2 (i.e., 
3) to have just one and not two independent random variables. An equivalent viewpoint is that each 
'wrapped' diagonal is periodic with period m and has m distinct random variables. Note that the 
diagonals are wrapped toroidally, and each such diagonal has N elements. 

Clearly if m = 1 these ensembles reduce to the previous cases, and as m — > oo they approach the 
full family of real symmetric matrices; in other words, the circulant or Toeplitz structure vanishes as 
m — > oo, but for any finite m there is additional structure. The goal of this paper is to determine the 
limiting spectral measures for these families and to quantify how the convergence to the semi-circle 
depends on m. We find an explicit closed form expression for the limiting spectral density of the 
m-block circulant family as a product of a Gaussian and a degree 2m — 2 polynomial. 

1.2. Results. Before stating our results, we must define the probability spaces where our ensemble 
lives and state the various types of convergence that we can prove. We provide full details for the 
m-block circulant matrices, as the related Toeplitz ensemble is similar. The following definitions 
and set-up are standard, but are included for completeness. We paraphrase from HMMSllJMPH with 
permission. 

Fix m and for each integer N let Q m ,N denote the set of m-block circulant matrices of dimension 
N. Define an equivalence relation ~ on {1, 2, ... , N} 2 . Say that (i, j) ~ (i', j') if and only if 
for all m-block circulant matrices, in other words, if 

• j — i = j' — i' (mod N) and i = i' (mod m), or 

• j — i = —(j' — i') (mod N) and i = j' (mod m). 

Consider the quotient {1,2,..., N} 2 {1,2, ... , N} 2 / ~. This induces an injection Rl 1 ' 2 '-.-^} 2 /- 
R^ 2 . The set Rl 1 ' 2 '-'^) 2 /- has the structure of a probability space with the product measure of 
p(x) dx with itself |{1, 2, . . . , N} 2 / ~ | times, where dx is Lebesgue measure. We define the prob- 
ability space (£l m> N, J- m ,N, Pm,7v) to be its image in R^ 2 = M N 2(R) under the injection, with the 
same distribution. 

To each A N E £l m ,N we attach a measure by placing a point mass of size 1 /n at each normalized 
eigenvalue Xi(A N ): 

N 



dx, 



(1.7) 



where S(x) is the standard Dirac delta function; see Footnoted] for an explanation of the normaliza- 
tion factor equaling \/N. We call [i m ,A N the normalized spectral measure associated with A^. 

Definition 1.2 (Normalized empirical spectral distribution). Let A N 6 Vt m ^ N have eigenvalues 
Ajv > ■ • • > Ai. The normalized empirical spectral distribution (the empirical distribution of 
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normalized eigenvalues) Fm N ^ is defined by 

= = V^'> . (1.8) 

As Fm N ^(x) = f* OQ fJ'm,A r f(t)dt, we see that Fm N ^ is the cumulative distribution function 

associated to the measure /i n ,A N - We are interested in the behavior of a typical Fm N ^ as we vary 
An in our ensembles il m ,N as N — > oo. 

Consider any probability space J7 m which has the f2 m) 7v as quotients. (The most obvious example 
is the independent product.) This paper build on a line of papers HHMl IMMSl I JMPH concerning var- 
ious Toeplitz ensembles which fix fi m to be the space of N-indexed strings of real numbers picked 
independently from p, with quotient maps to each Q m ,jv mapping a string to a matrix whose free 
parameters come from an initial segment of the right length. There is no need for the specificities 
of this construction, so we consider the general case. 

Definition 1.3 (Limiting spectral distribution). If as N — > oo we have Fm N ^ converges in some 
sense (for example, in probability or almost surely) to a distribution F m , then we say F m is the 
limiting spectral distribution of the ensemble. 

We investigate the symmetric m-block Toeplitz and circulant ensembles. We may view these as 
structurally weakened real symmetric Toeplitz and circulant ensembles. When m is 1 we regain the 
Toeplitz (circulant) structure, while if m = N we have the general real symmetric ensemble. If m 
is growing with the size of the matrix, we expect the eigenvalues to be distributed according to the 
semi-circle law, while for fixed m we expect to see new limiting spectral distributions. 

Following the notation of the previous subsection, for each integer N we let an d de- 
note the probability space of real symmetric m-block Toeplitz and circulant matrices of dimension 
N, respectively. We now state our main results. 

Theorem 1.4 (Limiting spectral measures of symmetric block Toeplitz and circulant ensembles). 

Let m 1 7V. 

(1) The characteristic function of the limiting spectral measure of the symmetric m-block cir- 
culant ensemble is 

K{t) = -e-^-e-^L^lJ-) = e- t2 / 2m M(m + l,2,-*y m ), (1.9) 
m \m J 

where L^__ l is a generalized Laguerre polynomial and M a confluent hypergeometric func- 
tion. The expression equals the spectral characteristic function for the m x m GUE. The 
limiting spectral density function (the Fourier transform of<p m ) is 



For any fixed m, the limiting spectral density is the product of a Gaussian and an even 
polynomial of degree 2m — 2, and has unbounded support. 
(2) If m tends to infinity with N (at any rate) then the limiting spectral distribution of the 
symmetric m-block circulant and Toeplitz ensembles, normalized by rescaling x to x/2, 
converge to the semi-circle distribution; without the renormalization, the convergence is to 
a semi-ellipse, with density /wig (see ( 11.31) ). 
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Figure 1 . Plots for fx, f 2 , f±, f$, fie and the semi-circle density. 

(3) As m — > oo, the limiting spectral measures f m of the m-block circulant ensemble converge 
uniformly and in L p for any p > 1 to /wig. with \f m (x) — fwig{ x )\ ^ m~ 2 ^ +e for any e > 0. 

(4) The empirical spectral measures of the m-block circulant and Toeplitz ensembles converge 
weakly and in probability to their corresponding limiting spectral measures, and we have 
almost sure convergence ifp is an even function. 

Figure Q] illustrates the convergence of the limiting measures to the semi-circle; numerical simu- 
lations (see Figures [2l [3] and HI) illustrate the rapidity of the convergence. We see that even for small 
m, in which case there are only mN /2 non-zero entries in the adjacency matrices (though these can 
be any of the iV 2 — N non-diagonal entries of the matrix), the limiting spectral measure is close to 
the semi-circle. This behavior is similar to what happens with (/-regular graphs, though in our case 
the convergence is faster and the support is unbounded for any finite m. 

Finally, the limiting eigenvalue density for m-block circulant matrices is the same as the eigen- 
value density of a certain Gaussian Hermitian ensemble. Specifically, we consider m x m Hermit- 
ian matrices with off-diagonal entries picked independently from a complex Gaussian with density 
function p(z) = ^e~' 2 ' , and diagonal entries picked independently from a real Gaussian of mean 
and variance 1. We provide a heuristic for why these densities are the same in S5.ll see also HZvH 
(especially Section 5.2) for a proof. 

Our results generalize to related ensembles. For example, the (wrapped) diagonals of our m- 
block circulant ensembles have the following structure (remember we assume m\N): 

)• (LID 

Note that we have a periodic repeating block of size m with m independent random variables; for 
brevity, we denote this structure by 

(di,d 2 ,...,d m ). (1.12) 
Similar arguments handle other related ensembles, such as the subfamily of period m-ciculant 
matrices in which some entries within the period are forced to be equal. Interesting comparisons 
are (d\, d 2 ) = (di, d 2 , di, d 2 ) versus (d\, di, d 2 , d 2 ) or (d±, d 2 , d 2 , d\). While it is a natural guess that 
the limiting spectral measure is determined solely by the frequency at which each letter appears, 
this is false. 
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Figure 2. (Left) Plot for f 2 and histogram of eigenvalues of 1000 symmetric period 
2-block circulant matrices of size 400 x 400. (Right) Plot for f 3 and histogram of 
eigenvalues of 1000 symmetric period 3-block circulant matrices of size 402 x 402. 




Figure 3 . (Left) Plot for / 4 and histogram of eigenvalues of 1000 symmetric period 
4-block circulant matrices of size 400 x 400. (Right) Plot for f 8 and histogram of 
eigenvalues of 1000 symmetric period 8-block circulant matrices of size 400 x 400. 




Figure 4. (Left) Plot for fi and histogram of eigenvalues of 1000 symmetric period 
1 -block circulant matrices of size 400 x 400. (Right) Plot for f 20 and histogram of 
eigenvalues of 1000 symmetric period 20-block circulant matrices of size 400 x 400. 



Theorem 1.5. Let V = (d^, d i2 , . . . , d im ) where each d, h e {d±, . . . , d u } and each di occurs 
exactly r, times in the pattern V, with r± + • ■ • + r v = m; equivalently, V is a permutation of 
(di, . . . , d±, d 2 , ■ ■ ■ , d 2 , ■ ■ ■ , d v , . . . , dy) with rj copies of di. Modify the N x N period m-block 
circulant matrices by replacing the pattern (di, d 2 , ■ ■ ■ , d m ) with V ( remember m\ N). Then for any 
V as N — oo the limiting spectral measure exists. The resulting measure does not depend solely on 
the frequencies of the letters in the pattern but also on their locations; in particular, while the fourth 
moments of the measures associated to {d±, d 2 , d±, d 2 } and {d±, di, d 2 , d 2 } are equal (interestingly, 
the fourth moment of any pattern only depends on the frequencies), the sixth moments differ. 
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We prove our main results using the method of moments. As the proof of Theorem II .51 is sim- 
ilar to that of Theorem 11.41 we just sketch the ideas and computations in Appendix |Bj For our 
ensembles, we first show that the average of the k th moments over our ensemble converge to the 
moments of a probability density. By studying the variance or fourth moment of the difference of 
the moments of the empirical spectral measures and the limits of the average moments, we obtain 
the various types of convergence by applications of Chebyshev's inequality and the Borel-Cantelli 
Lemma. These arguments are similar to previous works in the literature, and yield only the exis- 
tence of the limiting spectral measure. 

Unlike other works for related ensembles, however, we are able to obtain explicit closed form ex- 
pressions for the moments for the symmetric m-block circulant ensemble. This should be compared 
to the Toeplitz ensemble case, where previous studies could only relate these moments to volumes 
of Eulerian solids or solutions to systems of Diophantine equations. Similar to other ensembles, we 
show that the only contribution in the limit is when k — 21 and the indices are matched in pairs 
with opposite orientation. We may view this as a 2£-gon with vertices (ii, i 2 ), (12, h), ■ ■ ■ , (hi, h)- 
The first step is to note that when m = 1, similar to the circulant and palindromic Toeplitz ensem- 
bles, each matching contributes 1; as there are (2£ — 1)!! ways to match 2£ objects in pairs, and as 
(2£ — 1)!! is the 2£ th moment of the standard normal, this yields the Gaussian behavior. For general 
m, the key idea is to look at the dual picture. Instead of matching indices we match edges. In the 
limit as N — > oo, the only contribution occurs when the edges are matched in pairs with opposite 
orientation. Topologically, these are exactly the pairings which give orientable surfaces. If g is the 
genus of the associated surface, then the matching contributes m~ 29 . Harer and Zagier IIHarZall 
determined formulas for e g (£), the number of matchings that form these orientable surfaces. This 
yields the N — > oo limit of the average 2£ th moment is 



After some algebra, we express the characteristic function (which is the inverse Fourier transform; 
see Footnote© of the limiting spectral measure as a certain term in the convolution of the associated 
generating function of the e g 's and the normal distribution, which we can compute using Cauchy's 
residue theorem. Taking the Fourier transform (appropriately normalized) yields an explicit, closed 
form expression for the density. We note that the same formulas arise in investigations of the 
moments for Gaussian ensembles; see Section 1 .6 of IfFoll and HZvl (as well as the references therein) 
for additional comments and examples. 

The paper is organized as follows. In O we describe the method of proof and derive useful 
expansions for the moments in terms of quantities from algebraic topology. We use these in §[3] to 
determine the limiting spectral measures, and show convergence in §@] We conclude in §[5] with 
a description of future work and related results. Appendix lAl provides some needed estimates for 
proving the rate of convergence in Theorem 1 1.41 and we conclude in Appendix [B] with a discussion 
of the proof of Theorem 11.51 (see llXHl for complete details). 



In this section we investigate the moments of the associated spectral measures. We first describe 
the general framework of the convergence proofs and then derive useful expansions for the average 
moments for our ensemble for each N (Lemma [2T2T ). The average odd moments are easily seen to 




(1.13) 



2. Moments Preliminaries 
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vanish, and we find a useful expansion for the 2k th moment in Lemma |2~4l relating this moment to 
the number of pairings of the edges of a 2&;-gon giving rise to a genus g surface 

2.1. Markov's Method of Moments. For the eigenvalue density of a particular N x N symmetric 
m-block circulant matrix A, we use the redundant notation fi m ,A,N(x) dx (to emphasize the N 
dependence), setting 

N 

i=i x v 
To prove Theorem 1 1.4[ we must show 

(1) as iV — > oo a typical matrix has its spectral measure close to the system average; 

(2) these system averages converge to the claimed measures. 

The second claim follows easily from Markov's Method of Moments, which we now briefly 
describe. To each integer k > we define the random variable X k . mtN on Vt m by 



/j lAjN (x)dx := ^- ^2 5 ( x - ^j=- J dx. (2.1) 



X k , m , N {A) = / x k dF^(x); (2.2) 

J — oo 

note this is the k th moment of the measure n m ,A,N- 

Our main tool to understand the average over all A in our ensemble of the Fm N s is the 
Moment Convergence Theorem (see ItTall for example); while the analysis in HMMS1I was simplified 
by the fact that the convergence was to the standard normal, similar arguments (see also HJMPIO hold 
in our case as the growth rate of the moments of our limiting distribution implies that the moments 
uniquely determine a probability distribution. 

Theorem 2.1 (Moment Convergence Theorem). Let {Fn(x)} be a sequence of distribution func- 
tions such that the moments 

/CO 
x h dF N (x) (2.3) 
-oo 

exist for all k. Let {M^}^ be a sequence of moments that uniquely determine a probability 
distribution, and denote the cumulative distribution function by 7flim7v->oo -\ //,-.. v = M k then 
lim A r^ 00 F A r(x) = $(x). 

We will see that the average moments uniquely determine a measure, and will be left with proving 
that a typical matrix has a spectral measure close to the system average. The n th moment of A's 
measure, given by integrating x n against [i m ,A,N, is 



7 — 1 \ V / 7 — 1 



We define 

M n , m (N) := E(M n . m (A,N)), (2.5) 

and set 

M n . rn := lim M n . m (N) (2.6) 

7V->oo 

(we'll show later that the limit exists). By E(M n;m (A, N)), we mean the expected value of M n;m (A, N) 
for a random symmetric m-block circulant matrix A E £l m ,N- 
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2.2. Moment Expansion. We use a standard method to compute the moments. By the eigenvalue 
trace lemma, 

TV 

Tr(A n ) = (2-7) 
t=i 

so 

M n;m (A,N) = J J / ^Tr(A n ). (2.8) 

Expanding out Ti(A n ), 

M n . m (A,N) = N l /2+1 a ili2 a i2i3 ---a inil , (2.9) 

l<ti,...,i n <JV 

so by linearity of expectation, 

M n]m (N) = N } /2+1 E(a ili2 a i2i3 ---a inil ). (2.10) 

l<ii,...,i n <N 

Recall that we've defined the equivalence relation ~ on {1, 2, ... , iV} 2 by — {i',f) if and 

only if ciij = ai>ji for all real symmetric m-block circulant matrices. That is, — {i',f) if and 
only if 

• j — % = f — il (mod N) and % = il (mod m), or 

• j — i = — (f — i') (mod N) and i = j' (mod m). 

For each term in the sum in (12.101) . ~ induces an equivalence relation ~ on {(1, 2), (2, 3), . . . , 
(n, 1)} by its action on {(z'i, i 2 ), {12, iz), ■ ■ ■ , (i n , h)}- Let t?(~) denote the number of n-tuples with 
< ii, . . . , i n < N whose indices inherit ~ from ~. Say ~ splits up {(1, 2), (2, 3), . . . , (n, 1)} 
into equivalence classes with sizes c?i(~), . . . , <i;(~). Because the entries of our random matrices 
are independent identically distributed, 

E(a ili2 a i2i3 • • -a inil ) = m dl ^ym dl ^), (2.11) 

where the rrid are the moments of p. Thus, we may write 

M n]m {N) = N l /2+l V(~)m dli „) ■ ■ ■ m dlH . (2.12) 

As p has mean 0, m dl r^) ■ ■ ■ m dl (^) = unless all of the dj are greater than 1. So all the terms in 
the above sum vanish except for those coming from a relation ~ which matches at least in pairs. 

The r?(~) denotes the number of solutions modulo N the following system of Diophantine equa- 
tions: Whenever (s, s + 1) ~ (t, t + 1), 

• i s+ i — i s = z't+i — it (mod N) and i s = it (mod m), or 

• i 8+1 -i s = -(i t +i - k) (mod N) and i s = i t+1 (mod m). 

This system has at most 2 n ~ l N l+1 solutions, a bound we obtain by completely ignoring the 
(mod m) constraints (see also HMMSIO . Specifically, we pick one difference i s+ i — i s from each 
congruence class of ~ freely, and we are left with at most 2 choices for the remaining ones. Finally, 

we pick i\ freely, and this now determines all the i s — i\ + ^^(v+i — v). This method will not 

s'<s 

always produce a legitimate solution, even without the (mod m) constraints, but it suffices to give 
an upper bound on the number of solutions. 
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Z5 %<o 



Figure 5 . Diagram for a pairing arising in computing the 6 



moment. 



When n is odd, say n — 2k + 1, then / is at most k. Thus N nj 2 +i ' r ]( r ^) — N k+ 3 /2 ^ n 

-i N i+i < 

_t_2 n ~ l N k+l = -j=1 n ~ l = O n ( -4= ) . This implies the odd moments vanish in the limit, as 



When n is even, say n = 2k, then I is at most k. If / < k, then I < k — 1, and we have, similar 



to the above, j^tvH < J ^ T 2 n ~ l N l+1 < ^2 n ~ l N k = ±2 n ~ l = O n (i) . If Z = k, then the 



entries are exactly matched in pairs, that is, all the dj = 2. As p has variance 1 (i.e., m 2 = 1), the 
formula for the even moments, (12.121) . becomes 



We've changed notation slightly. The sum is now over pairings a on {(1, 2), (2, 3), . . . , (n, 1)}, 
which we may consider as functions (specifically, involutions with no fixed points) as well as equiv- 
alence relations. We have thus shown 

Lemma 2.2. For the ensemble of symmetric m-block circulant matrices, 



where the sum is over pairings a on {(1, 2), (2, 3), . . ., (n, 1)}. In particular, as N — > 00 the 
average odd moment is zero. 

2.3. Even Moments. We showed the odd moments go to zero like 1 /\/n as N — > 00; we now 



calculate the 2k th moments. From Lemma [2T2l the only terms which contribute in the limit are 
those in which the a^+i's are matched in pairs. We can think of the pairing as a pairing of the 
edges of a 2/c-gon with vertices 1,2, ... ,2k and edges (1, 2), (2, 3), ... , (2k, 1). The vertices are 
labeled i\, . . . , i 2 k and the edges are labeled a ili2 , . . . , a i2kil . See Figure[5] 

Note that this is dual to the diagrams for pairings that appear in HHMUMMSll . in which the ai 3 i s+1 
are represented as vertices. For more on such an identification and its application in determining 
moments for random matrix ensembles, see llFol (Section 1.6) and HZvL 

If a-isis+i an d a i t i t +i are paired, we have either 




(2.13) 




(2.14) 




(2.15) 
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Figure 6. Some possible orientations of paired edges for the 6-gon. 

• i s +i — i s = i t +i — it (mod N) and i s = i t (mod m), or 

• - i s = -(i t +i - h) (mod N) and z s = i t+1 (mod m). 

We think of these two cases as pairing (s, s+1) and (t, with the same or opposite orientation, 
respectively. For example, in Figure [6] the hexagon on the left has all edges paired in opposite 
orientation, and the one on the right has all but the red edges paired in opposite orientation. 

We now dramatically reduce the number of pairings we must consider by showing that the only 
pairings which contribute in the limit are those in which all edges are paired with opposite orienta- 
tion. Topologically, these are exactly the pairings which give orientable surfaces HHatllHarZal . This 
result and its proof is a minor modification of their analogs in the Toeplitz and palindromic Toeplitz 
cases HHCTllMMgllJMFl . 

Lemma 2.3. Consider a pairing o with orientations Ej. If any Ej is equal to 1, then the pairing 
contributes O^/n). 

Proof. The size of the contribution is equal to the number of solutions to the k equations 

i s+ i -i s = Ej(i a{s ) +1 - v (s) ) (mod N), (2.16) 

as well as some (mod m) equations, divided by N k+l . We temporarily ignore the (mod m) con- 
straints and bound the contribution from above by the number of solutions to the (mod N) equa- 
tions over iV fe+1 . Because the i s are restricted to the values 1, 2, . . . , iV, we can consider them as 
elements of Z/NZ, and we now notate the (mod N) congruences with equality. 

The pairing puts the numbers 1,2, ... ,2k into k equivalence classes of size two; arbitrarily order 
the equivalence classes and pick an element from each to call Sj, naming the other element tj = 

Our Z/NZ equations now look like 

i 8j+1 - i Sj = ej(i tj+1 - i tj ) mod N. (2.17) 

Defining 

Uj '■= Hj+i ~ Hji 
our equations now look like Xj = EjUj. Thus 

= J2 is + i ~ is = X^'+X^ 

s=i i=i i=i 



(2.18) 

k 

i=i 
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If any one of the Ej = 1, this gives a nontrivial relation among the yj, and we lose a degree of 
freedom. We may choose k — 1 of the jjj freely (in Z/JVZ), and we are left with 1 or possibly 2 
choices for the remaining y,j (depending on the parity of N). The Xj are now determined as well, 
so z s+ i — i s is now determined for every s. If we choose i\ freely, this now determines all the 
i s = iy + y^(j s , +1 - v). Thus, we have at most A^ -1 • 2 • N = 2N k solutions to (12.161) . So the 

s'<s 

contribution from a pairing with a positive sign is at most Ok(2 Nk /N k+1 ) = O^/n). (The reason 
for the big-Oh constant depending on k is that if some of the different pairs have the same value, 
we might not have k copies of the second moment but instead maybe four second moments and two 
eighth moments; however, the contribution is trivially bounded by maxi<g<fe(l + m 2 ^) fc , where m 2 ^ 
is the 2£ th moment of p.) □ 

Thus we have 

M 2k , m {N) = J2^)N~ {h+1) + O k , (2.19) 

where in (a) denotes the number of solutions to 

ij+i - ij = -(iaj+i - iaj) mod N (2.20) 

and 

ij = i*(j)+i, ij+i = mod m (2.21) 

(the second (mod m) constraint is redundant). We discuss how to evaluate this moment in closed 
form, culminating in Lemma l2~4l 

We now consider a given pairing as a topological identification (see HHatll for an exposition of the 
standard theory); this is the crux of our argument. Specifically, consider a 2A;-gon with the interior 
filled in (homeomorphic to the disk), and identify the paired edges with opposite orientation. Under 
the identification, some vertices are identified; let v denote the number of vertices in the quotient. 

Consider the (Z/AZ)-submodule A of (Z/iVZ) 2fc in which the (mod N) constraints hold. We 
have A is isomorphic to (Z/iVZ) fe+1 . Specifically, we may freely choose the value of exactly half 
of the differences i s+ \ — i s , and then the rest are determined. Because all the pairings are opposite 
orientation, these "differences" sum to zero, so they are actually realizable as differences. Now 
choose i\ freely, and the rest of the i s = i\ + ^^(v+i — i s >) are determined. 

s'<s 

Let A denote the quotient of A in which everything is reduced modulo m, and consider the 
(Z/mZ)-submodule B C A in which the modulo m constraints hold. By (12.211) . we can see that 
the labels at two vertices of our 2A;-gon are forced to be congruent (mod m) if and only if the 
vertices are identified in the quotient, and these are all the (mod m) constraints. In other words, 
B is isomorphic to (Z/mZ) 1 '. An element of A for which the (mod m) constraints also hold is 
exactly one in the preimage of B. We have m v choices for an element in B, and there are ( N /m) k+1 
ways to lift such an element to an element of A in its fiber. Thus, the equations have a total of 
m v ( N /m) k+1 = m ~( k+1 ~ v ' ) N k+1 , so the pairing has a contribution of m~^ k+1 ~ v \ 

Let X be the 2-dimensional cell complex described by the pairing a of the edges of the 2k- 
gon. Because all edges were paired in the reverse direction, X is an orientable surface. After 
identifications, the complex we've described has 1 face, k edges, and, say, v vertices. If we denote 
by g the genus of the surface, we obtain two expressions for the Euler characteristic of X. By the 
standard (homological) definition of Euler characteristic, we have xPO = 1 — k + v. On the other 
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hand, for a genus g surface X, x(X) = 2 — 2g HHatL Equating and rearranging, 

2g = k + l- v. (2.22) 
Thus the pairing a contributes m~ 29 , and we have shown 
Lemma 2.4. For the ensemble of symmetric m-block circulant matrices, 

M 2k , m (N) = J2 £ 9 ( k ) m ~ 29 + °k (^) , (2-23) 

9 

where s g (k) denote the number of pairings of the edges of a 2k-gon which give rise to a genus g 
surface. 

3. DETERMINING THE LIMITING SPECTRAL MEASURES 

We prove parts (1) and (2) of Theorem 1 1.41 Specifically, we derive the density formula for the 
limiting spectral density of symmetric m-block circulant matrices. We show that, if m grows at any 
rate with N, then the limiting spectral density is the semi-circle for both the symmetric m-block 
circulant and Toeplitz ensembles. 

3.1. The Limiting Spectral Measure of the Symmetric m-Block Circulant Ensemble. 

Proof of Theorem [L4\ l ). By deriving an explicit formula, we show that the limiting spectral den- 
sity function f m of the real symmetric m-block circulant ensemble is equal to the spectral density 
function of the m x m GUE. 

From Lemma [2~4l the N — > oo limit of the average 2k th moment equals 

LVaJ 

M 2k]m = J2 £ M m ~ 29 , (3-D 

9=0 

with e g (k) the number of pairings of the edges of a 2A;-gon giving rise to a genus g surface. Harer 
and Zagier HHarZall give formulas for the £ g (k). They prove 

eJk) = ^ -x ( coefficient of x 2g in ( — r^-rx) | (3.2) 

gK 1 (k + l)\(k - 2g)\ \ \tmh{*/2)) ) 

and 

LVaJ 

J2^g(k)r k+1 ~ 29 = (2k-l)\\c(k,r), (3.3) 

9=0 



where 



Thus, we may write 



1 + 2J2 c(k,r)x k+1 = I til) . (3.4) 

k=0 ^ X ' 



M 2k , m = m- ( - k+1 \2k-l)\\c(k,m). (3.5) 
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We construct the characteristic functiorOof the limiting spectral distribution. Let X m be a random 
variable with density f m . Then (remembering the odd moments vanish) 



itx , ^ (ttYM e , m 



t (t) = E[e*«»] = 

m 2k M 2k ., 



oo 

^ 7 ^m-"= +1, (2<:-l)!!c(*,m)(-i 2 )*. (3.6) 



o (»)' 

In order to obtain a closed form expression, we rewrite the characteristic function as 

1 00 1 ( ' -t 2 \ k 

*.«) = ^E*-) M (^j, (3-7) 

fc=0 v 7 

using (2A; — 1)!! = |St. The reason for this is that we can interpret the above as a certain coefficient 
in the convolution of two known generating functions, which can be isolated by a contour integral. 
Specifically, consider the two functions 

i / / 1 I \ m \ 00 00 fc 

F(y):= ^ ((rrf - 1 = E c ^-)^ and ^) : = ey = E|r- < 3 - 8) 

Note that <p m {t) is the function whose power series is the sum of the products of the k th coefficients 
of G(—y 2 /2m) (which is related to the exponential distribution) and F(y) (which is related to the 
generating function of the s g (k)). Thus, we may use a multiplicative convolution to find a formula 
for the sum. By Cauchy's residue theorem, integrating F(z~ 1 )G(— t2z /2m)z~ 1 over the circle of 
radius 2 yields 

since the constant term in the expansion of F(z~ 1 )G(— t2z /2m) is exactly the sum of the products of 
coefficients for which the powers of y in F(y) and G{y) are the sameQ We are integrating along the 



3 The characteristic function is (f) m {t) = IL[e itXm ] = f m {x)e ltx dx. This is the inverse Fourier transform of 

fm- 

4 A11 functions are meromorphic in the region with finitely many poles; thus the contour integral yields the sum of 
the residues. See for example MSS2I . 
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circle of radius 2 instead of the unit circle to have the pole inside the circle and not on it. Thus 



2irim J\ z \ =2 2z~ x \ \ 1 — z~ 
1 f ffz + 1 



Anim J ]z]=2 \\z-l 



1 e- t2z ' 2m dz 



t 2 /2r, 



1 + -=- -1 ) e-^-^dz 



A-nim J\ z \ =2 \ \ z-l 

e -^m , - /m w 2 Vol f-e 



A-nim /ui-9 \ I / V z — 1 / ^— ' s! \ 2m 



3 t2(z - 1] / 2m dz. (3.10) 



47rim yj z(=2 

By Cauchy's Residue Theorem the second integral vanishes and the only surviving terms in the 
first integral are when I — s = 1, whose coefficient is the residue. Thus 

= 1 e- a/2 ™£ (?) TT^ (-)"= 1 ^-e, ('%.) , (3.1D 
m \ £ j [l — 1)1 \ m J m 

which equals the spectral density function of the m x m GUE (see ULedll ). 

As the density and the characteristic function are a Fourier transform pair, each can be recovered 
from the other through either the Fourier or the inverse Fourier transform (see for example USSll 
ISS2I0 . Since the characteristic function is given by 



f oo 

<j> m (t) = E[e ltx ™] = I e ltx f m {x)dx (3.12) 

J — oo 

(where X m is a random variable with density f m ), the density is regained by the relation 

I r°° 

f m (x) = Kix) = — / e~ ltx <p m {t)dt. (3.13) 

Taking the Fourier transform of the characteristic function (j> m (t), and interchanging the sum and 
the integral, we get 



- 4 ECMui ( - mr ' '-■ <3J4) 
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Completing the square in the integrand of I m , we obtain 

I m = e~ mx2/2 j°° t^ l ~ l) exp (~ + i^x^J ) dt. (3.15) 

Changing variables by y = -^=t + iy/mx, dy = -j^dt, we find I m equals 

f oo 



/oo 
(y - i^ix) 2{l ~ l) (V^) 2{1 ~ 1] e'^V^dy 
-oo 

2 J2 ( 2(l " (-^) 2{l - iyS r y s e~^dy. (3.16) 

_n V s / J -ex, 



-mx 2 /2 I— 



2(1-1) 

_1 

"m" 2 

s=0 

The integral above is the s th moment of the Gaussian, and is \/2n(s — 1)!! for even s and for odd 
s. Since the odd s terms vanish, we replace the variable s with 2s and sum over < s < (I — 1). 
We find 

I m = v^e-^m^ £ ( 2{l ~ l A (-mx 2 ) 1 - 1 - 3 (2s - 1)!!. 

(3.17) 

Substituting this expression for I m into (13.141) and making the change of variables r — I — 1 — s, 
we find that the density is 

/ml ' VSS ^ (2r)! \r + s + l) (r + ,)!,! ( 2 J J {mX > ' (3 ' 18) 

This completes the proof of Theorem 11.41 1 ) . □ 

3.2. The m — > oo Limit and the Semi-Circle. Before proving Theorem II .41 2), we first derive 
expressions for the limits of the average moments of the symmetric m-block Toeplitz ensemble. 
We sketch the argument. Though the analysis is similar to its circulant cousin, it presents more 
difficult combinatorics. Because diagonals do not "wrap around", certain diagonals are better to be 
on than others. Consequently, the Diophantine obstructions of HHMH are present. The problems are 
the matchings with "crossings", or, in topological language, those matchings which give rise to tori 
with genus g > 1 as opposed to spheres with g = 0. For a detailed analysis of the Diophantine 
obstructions and how the added circulant structure fixes them, see HHMH and HMMSH . Fortunately, 
it is easy to show that the contributions to the 2k th moment of the symmetric m-block Toeplitz 
distribution from the non-crossing (i.e, the spherical matchings or, in the language of HBanBoL the 
Catalan words) are unhindered by Diophantine obstructions and thus contribute fully. The number 
of these matchings is C fe , which is the /c* Catalan number ^j-j- ( 2 fc fc ) as well as the 2k^ moment of 
the Wigner density 

/wig(x) = iW 1_ (i) 2 if W^ 2 (3.19) 
[ otherwise. 

Note that with this normalization have a semi-ellipse and not a semi-circle; to obtain the semi- 
circle, we normalize the eigenvalues by 2\/N and not \/N. As the other matchings contribute zero 
in the limit, we obtain convergence to the Wigner semi-circle as m — > oo. We now prove the above 
assertions. 
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Lemma 3.1. The limit of the average of the 2k moment of the symmetric m-block Toeplitz ensem- 
ble equals 

|k/3j 

M 2k , m = C k + J2 d ( k ^)m- 29 , (3.20) 

9=1 

where C k is the k th Catalan number and d(k, g) G [0, 1] are constants corresponding to the total 
contributions from the genus g pairings for the 2k th moment. 

Proof. For the symmetric m-block Toeplitz ensemble, the analysis in §|2] applies almost exactly. In 
the condition for = a^ji, equality replaces congruence modulo N. 

• j — % — f — i' and % = %' (mod m), or 

• j — i = — (j' — i') and i = j' (mod m). 

These constraints are more restrictive, so we again obtain 2 n ~ l N l+1 as an upper bound on the 
number of solutions. Following the previous argument, the odd moments are M 2k+ i- m (N) = 
OkiyVw), and the even moments are 

M 2fc;m (A0 = E lW + °* (jf) > ( 3 - 21 ) 

where rj(a) is the number of solutions to the Diophantine equations arising from the pairings a on 
{(1, 2), (2, 3), ... , (2k, 1)} of the indices. Thus the odd moments vanish in the limit. Moreover, the 
only matchings that contribute are the ones with negative signs. To see this fact, one can follow the 
proof of Lemma [231 except working in Z instead of Z/7VZ. 

While it is known that most matchings for the real symmetric Toeplitz ensemble do not contribute 
fully, a general expression for the size of the contributions is unknown, though there are expressions 
for these in terms of volumes of Eulerian solids (see HBDJIO or obstructions to Diophantine equa- 
tions (see HHMIO . These expressions imply that each matching contributes at most 1. We introduce 
constants to denote their contribution (this corresponds to the m = 1 case). This allows us to handle 
the real symmetric m-block Toeplitz ensemble, and (arguing as in the proof of Lemma 12.41) . write 
the limit of the average of the 2k th moments as 

M 2k , m = J2 d ( k '9)m~ 2g . (3.22) 

9=0 

Here d(k, g) is the constant corresponding to the contributions of the genus g matchings. All that 
is left is to show that d(k, 0), the contributions from the non-crossing or spherical matchings, is the 
Catalan number C k . 

We know that the number of non-crossing matchings of 2k objects into k pairs is the Catalan 
number C k . This is well-known in the literature. Alternatively, we know the number of non- 
crossing matchings are e (k), as these are the ones that give the genus sphere. The claim follows 

immediately from (13.21 ) by taking the constant term (as g — 0) and noting tanh(|) = f — fj H • 

We are thus reduced to proving that, even with the mod m periodicity, each of these pairings still 
contributes 1. 

One way of doing this is by induction on matchings. Consider a non-crossing configuration of 
contributing matchings for the 2k th moment. Consider an arbitrary matching in the configuration, 
and denote the matching by a\. The matching corresponds to an equation i s — i s+1 = i t+1 — i t . 
If the matching is adjacent, meaning s = t + 1, then z f+1 is free and i t = i t+2 , and there is no 
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"penalty" (i.e., a decrease in the contribution) from the (mod m) condition. We call this having the 
ends of a matching "tied" (note that adjacent matchings always tie their ends). Otherwise, note that 
since we are looking at even moments, there are an even number of indices. Thus, to either side of 
the matching ol\ there can only be an even number of indices matched between themselves, since 
otherwise some matching would be crossing over a±. Thus, to either side, we are reduced to the 
non-crossing configurations for a lesser moment. By induction, these two sub-configurations are 
tied, and then trivially tie with our initial matched pair. As at each step there were no obstructions 
on the indices, this matching contributes fully, completing the proof. □ 

Our claims about convergence to semi-circular behavior now follow immediately. 

Proof of Theorem \L4\ 2 ). It is trivial to show that the symmetric m-block circulant ensemble has its 
limiting spectral distribution converge to the semi-ellipse as m — > oo because we have an explicit 
formula for its moments. From Lemma (12.41). we see that 



lim M 2k . m (N) = lim V = £ °( k )> (3-23) 

9< k h 

which in the proof of Lemma I3T1 we saw equals the Catalan number C k . 

We now turn to the symmetric m-block Toeplitz case. The proof proceeds similarly. From 
Lemma I3TI we have 



lim M 2k , m = lim \C k + V = Ck , (3.24) 

\ 9< k h / 

completing the proof. □ 



4. Convergence of the Limiting Spectral Measures 

We investigate several types of convergence. 

(1) (Almost sure convergence) For each k, X k . miN — > X kjTn almost surely if 

P m ({A e n rn : X k . rniN (A) -> X kjm (A) as iV — >■ oo}) = 1; (4.1) 

(2) (Convergence in probability) For each k, X k . mjN — > X k>m in probability if for all e > 0, 

lim F m (\X k . m>N (A) - X Km {A)\ > e) = 0; (4.2) 

N— >oo 

(3) (Weak convergence) For each k, X k . rnjN — > X k:Tri weakly if 

F m (X k . mtN {A) < x) -> ¥{X k , m {A) < x) (4.3) 

as N — > oo for all x at which Fx k m (x) := F(X k)m (A) < x) is continuous. 

Alternate notations are to say either with probability 1 or strongly for almost sure convergence and 
in distribution for weak convergence; both almost sure convergence and convergence in probability 
imply weak convergence. For our purposes we take X k>m as the random variable which is identically 
Mk,m, the limit of the average m th moment (i.e., limA^oo M kjm . N ), which we show below exist and 
uniquely determine a probability distribution for our ensembles. 

We have proved the first two parts of Theorem 1 1.41 which tells us that the limiting spectral mea- 
sures exist and giving us, for the symmetric m-block circulant ensemble, a closed form expression 
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for the density. We now prove the rest of the theorem, and determine the various types of conver- 
gence we have. We first prove the claimed uniform convergence of part (3), and then discuss the 
weak, in probability, and almost sure convergence of part (4). 

We use characteristic functions and Fourier analysis to show uniform (and thus pointwise) con- 
vergence of the limiting spectral distribution of the symmetric m-block circulant ensemble to the 
semi-ellipse distribution (remember it is an semi-ellipse and not a semi-circle due to our normal- 
ization). We note that this implies L p convergence for every p. The proof follows by showing the 
characteristic functions are close, and then the Fourier transform gives the densities are close. 

Proof of Theorem \L4\ 3 ). The density f m is the Fourier transform of 4> m (equivalently, m is the 
characteristic function associated to the density f m , where we have to be slightly careful to keep 
track of the normalization of the Fourier transform; see (13.121 )); similarly the Wigner distribution 
fwi g (x) is the Fourier transform of 0, where the Wigner distribution (a semi-ellipse in our case due 
to our normalizations) is 



jh\P--(W ifM<2 fAA , 
/wig(z) = \ 11 v K2J (4.4) 
otherwise. 

As our densities are nice, we may use the Fourier inversion formula to evaluate the difference. 
We find for any e > that 



|? m (x) - <f>(x) 



1 f°° 

— / (4> m (t)-(f>(t))e- ltx dt 

^ J-oo 

/oo 
\<f> m {t)-(f>{t)\dt 
-oo 

2/ 



< m" 2/9+e , (4.5) 

where the bound for this integral is proved in Lemma IA.1I and follows from standard properties 
of Laguerre polynomials and Bessel functions. Thus, as m — > oo, f m (x) = 4> m (x) converges to 
fwig(x) = 4>{x) for all x G 1. As the bound on the difference depends only on m and not on x, the 
convergence is uniform. 

We now show LP convergence. We have L°° convergence because it is equivalent to a.e. uni- 
form convergence. For 1 < p < oo, we automatically have L p convergence as we have both L 1 
convergence and the L°° norm is bounded. □ 

Proof of Theorem [L4\ 4 1 The proofs of these statements follow almost immediately from the argu- 
ments in HHM1 IMMSl IJMPII , as those proofs relied on degree of freedom arguments. The additional 
structure imposed by the (mod m) relations does not substantially affect those proofs (as can seen 
in the generalizations of the arguments from HHMM to HMMSI to HJMPII ). □ 

5. Future Research 

We discuss some natural, additional questions which we hope to study in future work. 

5.1. Representation Theory. The N x N m-block circulant matrices form a semisimple algebra 
over R. This algebra may be decomposed into N simple subalgebras of dimension m 2 , all but one 
or two of which are isomorphic to M m (C). One can show that, up to first order, this decomposition 
sends our measure on symmetric m-block circulant matrices to the m x m Gaussian Unitary Ensem- 
ble. One may then give a more algebraic proof of our results and circumvent the combinatorics of 
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Figure 7 . Density of nonzero spacings of the 10 central eigenvalues of 100 1024 x 
1024 symmetric circulant matrices, with independent entries picked i.i.d.r.v. from a 
Gaussian, normalized to have mean spacing 1. Compared to exponential and GOE 
densities. 

pairings; combining the two proofs gives a new proof of the results of HHarZaL This approach will 
appear in a more general setting in an upcoming paper of Kopp. The general result may be regarded 
as a central limit theorem for Artin-Wedderburn decomposition of finite-dimensional semisimple 
algebras 

5.2. Spacings. Another interesting topic to explore is the normalized spacings between adjacent 
eigenvalues. For many years, one of the biggest conjectures in random matrix theory was that 
if the entries of a full, N x N real symmetric matrix were chosen from a nice density p (say 
mean 0, variance 1, and finite higher moments), then as N — > oo the spacing between normalized 
eigenvalues converges to the scaling limit of the GOE, the Gaussian Orthogonal Ensemble (these 
matrices have entries chosen from Gaussians, with different variances depending on whether or not 
the element is on the main diagonal or not). After resisting attacks for decades, this conjecture was 
finally proved; see the work of Erdos, Ramirez, Schlein, and Yau HERSYl IESYH and Tao and Vu 



While this universality of behavior for differences seems to hold, not just for these full ensembles, 
but also for thin ensembles such as rf-regular graphs (see the numerical observations of Jakobson, 
(S. D.) Miller, Rivin and Rudnick HJ MRRfl ), we clearly do not expect to see GOE behavior for all 
thin families. A simple counterexample are diagonal matrices; as N — > oo the density of normalized 
eigenvalues will be whatever density the entries are drawn from, and the spacings between normal- 
ized eigenvalues will converge to the exponential. We also see this exponential behavior in other 
ensembles. It has numerically been observed in various Toeplitz ensembles (see HHMllMMSlI ). 

For the ensemble of symmetric circulant matrices, we cannot have strictly exponential behavior 
because all but 1 or 2 (depending on the parity of N /m) of the eigenvalues occur with multiplicity 
two. This can be seen from the explicit formula for the eigenvalues of a circulant matrix. Thus, the 
limiting spacing density has a point of mass | at 0. Nonetheless, the nonzero spacings appear to be 
distributed exponentially; see Figure[7J 

Similarly, for a symmetric m-block circulant matrix, all but N—m or N—m—1 of the eigenvalues 
occur with multiplicity two. The nonzero spacings appear to have the same exponential distribution 
(see Figured]). This is somewhat surprising, given that the eigenvalue density varies with m and 
converges to the semi-circle asm — > oo. While we see new eigenvalue densities for m constant, 
numerics suggest that we'll see new spacing densities for N /m constant. 



irrvniTvzi . 
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Figure 8. Density of nonzero spacings of the 10 central eigenvalues of 100 
1024 x 1024 symmetric m-block circulant matrices, with independent entries 
picked i.i.d.r.v. from a Gaussian, normalized to have mean spacing 1, with m = 
2, 16, 128, 256, 512, 1024, respectively. Compared to exponential and GOE densi- 
ties. 

However, for symmetric m-block Toeplitz matrices, we see different behavior (see Figure©. The 
spacings look exponentially distributed for m = 1 and appear to converge to the GOE distribution 
as we increase m. In the Toeplitz case, but not in the circulant, we see the spacings behaving as the 
spectral densities do. 

The representation theoretic approach will be used to solve the spacings problem for symmetric 
m-block circulant matrices in an upcoming paper of Kopp. The spacing problem for block Toeplitz 
matrices will require some new innovation. 

Appendix A. Pointwise Convergence as m -> oo 
This appendix by Gene Kopp, Steven J. Miller and Frederick Strauch^ 
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Figure 9. Density of spacings of the 10 central eigenvalues of 100 1024 x 1024 
symmetric m-block Toeplitz matrices, with independent entries picked i.i.d.r.v. from 
a Gaussian, normalized to have mean spacing 1, with m = 1,2,4, 16, 128, 1024, 
respectively. Compared to exponential and GOE densities. 



The characteristic function for the spectral measures of the period m-block circulant matrices is 



= 1 e"^ J2 



m 



1 



I ) (£-1)1 



-t l 



m 



(A.l) 



which solves the differential equation 

*C(t) + 30' m (t)+t(4 



<l>m{t) 







(A.2) 



with initial condition m (O) = 1; letting m — > oo gives t(f>"(t) + 30' (i) + 4t(f)(t) = 0, with initial 
condition 0(0) = 1. The solution to the finite m equation is a Laguerre polynomial, and the m = oo 
limit is J i( 2 *)/t with J\ the Bessel function of order 1. 



To see this, recall that the generalized Laguerre polynomial (see HASH ) has the explicit represen- 
tation 



■T) 



(A.3) 
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To compare (IA.1I) with (IA.3I) . we first shift the summation index by one (£ h-» t + 1) to find 

m-l / \ -i / j.2 \ ^ 



Using the identity 



m \ I m 

+ l) ~ U- 1 



(A.5) 



we see that n — m — 1, a = 1, and thus the characteristic function can be written in terms of the 
Laguerre polynomial: 

Kit) = -e~ t2 ^L^U t2 /m), (A.6) 
m 

or equivalently in terms of the confluent hypergeometric function 

= e- t2/2m M(m + l,2,-* 2 / m ). (A.7) 

From 13.2.2 of [ASJ we have limm^oo 4> m (t) = (f>(t); however, we need some control on the rate 
of convergence. 

Lemma A.l. Let r > 1 /3 and ft — |(1 — r). For all m and all t we have 

yt-# « r (A.8) 
It 11 + m /4 exp(— */2m) otherwise, 

where the implied constant is independent ofm but may depend on r. This implies 

\(p m (t)-(j>(t)\dt < m-V (A.9) 

Letting e > and taking r = ~ + 3e implies the integral is 0(m~ 2 ^ +t ). 

Proof. We first consider small t: \t\ < vnP with /3 = |(1 — r). Using 13.3.7 of HASH with a = m + 1, 
6 = 2 and z = — <2 /m to bound the confluent hypergeometric function M, we find 

<M*) = e- t2 / 2m M(m + l,2,-*7 m ) = + A n {2m)- n {-l) n t n ~ l J n+1 (2t), (A. 10) 

n=l 

where A = 1, A x = 0, A 2 = 1 and A n+i = A n -i + ^A n _ 2 for n > 2. 

For any r > l /z and m sufficiently large we have A n < m rn (we can't do better than r > l Jz 
as A 3 = |m). This follows by induction. It is clear for n < 2, and for larger n we have by the 
inductive assumption that 

A n+1 = A n _i + -^A n _ 2 < m r(n " 1} + m ■ m r(n " 2) = m r(n+1) • (m~ 2r + m 1 " 3 "); (A.ll) 
n + 1 

as r > !/3 the above is less than m r ^ n+1 ^ for m large. If we desire a bound to hold for all m, we 
instead use A n < c r m rn for c r sufficiently large. Substituting this bound for A n into (IA. 101) . noting 
Jim/t = <f>(t) and using | J n (x)\ < 1 (see 9.1.60 of flASl) yields, for \t\ < m 1 ~ r , 

oo / , \ n— 1 

< ^bEU^ <rm-^). (A.12) 

n=l ^ ' 
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We now turn to t large: \t\ > mP . Using 

\<p m (t) - <p(t)\ < \<f> m (t)\ + \<Kt)\ (A.13) 

to trivially bound the difference, the claim follows the decay of the Bessel and Laguerre functions. 
Specifically, (see 8.451(1) of llGRl ) we have J x {x) < x~ 1/2 and thus 

0(t) = < t - z '\ {AAA) 
t 

For <pm(t), we use 8.978(3) of E3D, which states 

L^(x) = Tr-V'e^ar-^-^n^-^cos (2^-^ - ^) + O (r^ 2 " 3 / 4 ) , (A.15) 

so long as Im(a) = and x > 0. Letting x = * 2 /m with |t| > h log 1 / 2 m, a = 1 and n = m — 1 we 
find 

m (t) = m^e-^L^/m) 



t2 /2m/f2/ \— 3/4 1/4 1 — 1/4 

e 7 (* /m) 7 m 7 + m ' 



< r 3 / 2 + m- 5/4 e- t2/2m . (A. 16) 

All that remains is to prove the claimed bound for \(f) m (t) — 4>(t) \ dt. The contribution from 
\t\ < mP is easily seen to be O r ( m/3 /m 1 - r ) = O r (m~ (1 ~ r)//3 ) with our choice of 0. For \t\ > mP , we 
have a contribution bounded by 



2 



00 . . POO 



t~ 3/2 + m" 5/4 e" t2/2m ) < m- fi/2 + m^ 3/4 / —== exp(-*7 2m )^ 

< m- (1_r)/3 + m" 3/4 , (A. 17) 



as the last integral is that of a Gaussian with mean zero and variance m and hence is 1 . (We chose 
f3 — |(1 — r) to equalize the bounds for the two integrals.) □ 



Appendix B. Generalized m-BLOCK Circulant Matrices 
This appendix by Steven J. Miller and Wentao Xion^ 

As the proofs are similar to the proof for m-block circulant matrices, we just highlight the differ- 
ences. The trace expansion from before holds, as do the arguments that the odd moments vanish. 

We first explore the modulo condition to compute some low moments, and show that the differ- 
ence in the modulo condition between the m-block circulant matrices and the generalized m-block 
circulant matrices leads to different values for moments, and hence limiting spectral distributions. 
Thus the limiting spectral distribution depends on the frequency of each element, as well as the way 
the elements are arranged, in an m-pattern. 
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Figure 10. The four zones for m-block circulant matrices. 



B.l. Zone- wise Locations and Pairing Conditions. Since we have restricted the computation of 
moments to even moments, and have shown that the only configurations that contribute to the 2k th 
moment are those in which the 2k matrix entries are matched in k pairs in opposite orientation, we 
are ready to compute the moments explicitly. We start by calculating the 2 nd moment, which by 
(12.101) is 2~2i<ij<N a ij a ji- As long as the matrix is symmetric, = aji and the 2 nd moment is 
1. We now describe the conditions for two entries a ia i 3+1 , a itit+1 to be paired, denoted as a isis+1 = 
a ith+i < ^ = ^ > ( s ? s + 1) ~ {t,t + l), which we need to consider in detail for the computation of higher 
moments. To facilitate the practice of checking pairing conditions, we divide an iV x N symmetric 
m-block circulant matrix into 4 zones (see Figure [TOl). and then reduce an entry a isis+1 in the matrix 
to its "basic form". Write it = mi]t + ee, where r)£ £ {1, 2, ... , ^} and q £ {0, 1, . . . , m — 1}, we 
have 

(1) < i s+1 - i s < f - 1 a lsis+1 £ zone 1 and a isls+1 = a esMvs+1 _ Vs)+es+1 ; 

(2) f < i s+1 - i s < N - 1 a lgls+1 £ zone 2 and a isis+1 = a es+iMll3+ ^ V3+l)+£a ; 

(3) f < i s - is+i <N-1^ a lals+1 £ zone 3 and a iats+1 = a 6a)TO( ^ +1+ ^_ IJ>)+6j+1 ; 

(4) < i s - i s+1 < f - 1 a isis+1 £ zone 4 and a isis+1 = a 6s+lMr]s _ Vs+l)+es . 

In short, (i s +i — i s ) determines which diagonal a is i s+1 is on. If a isi3+1 is in zone 1 or 3 (Area I), e s 
determines the slot of a isis+1 in an m-pattern; if a isis+1 is in zone 2 ox A (Area II), e s+ \ determines 
the slot of a isis+1 in an m-pattern. 



Recall the two basic pairing conditions, the diagonal condition that we have explored before, and 
the modulo condition, for which we will define an equivalence relation 1Z. For a real symmetric 
m-block circulant matrix following a generalized m-pattern and any two entries a is is+1 , a itit+1 in the 
matrix, suppose that i s and i t +i are the indices that determine the slot of the respective entries, then 
i s lZi t +i if and only if di s i s+1 , aj t i t+1 are in certain slots in an m-pattern such that these two entries can 
be equal. For example, for the {a, b} pattern, i s lZi t +\ < ^=^ i s = h+i (mod 2); for the {a, a, b, b} 
pattern, i s 1Zi t+ i ^=>- mod(z s ,4), mod(z m ,4) £ {1,2} or mod(z s ,4), mod(z 4+ i,4) £ 
{3,0}. 

We now formally define the two pairing conditions. 

(1) (diagonal condition) i s — i s+1 = — (i t — i t +i) (mod N). 

(2) (modulo condition) i s TZi t+ i or i s+1 lZi t , depending on which zone(s) a>i s i s+1 ,a>i t i t+1 are lo- 
cated in. 
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Since the diagonal condition implies a Diophantine equation for each of the k pairs of matrix entries, 
we only need to choose k + 1 out of 2k i/s, and the remaining i/s are determined. This shows 
that, trivially, the number of non-trivial configurations is bounded above by N k+1 . In addition, the 
diagonal condition always ensure that a isia+1 and a itit+1 are located in different areas. For instance, 
if a-isis+i e zone 1 and is ~ is+i = -(it - it+i), then a isis+1 e zone 4; if a isia+l e zone 1 and 
i s — i s+ i = — (i t — it+i) — N, then a isia+1 E zone 2, etc. Thus, if i s determines the slot for 
a i a i s +i m an m pattern, then i t+ i determines for a» t i t+1 ; if i s +i determines the slot for ai s i s+1 , then i t 
determines for a itit+1 , and vice versa. 

Considering the "basic" form of the entries, the two conditions above are equivalent to 

(1) (diagonal condition) (m?7 s + e s )-(m?7 s+1 + e s+ i) = - (mr) t + e t ) + (mr] t+1 + e t+1 ) (mod N) 

m(r] s - r] s+1 + 7] t - 7] t +i) + (e s - e s+1 + e t - e t+1 ) = or ±N. 

(2) (modulo condition) e s 1Ze t+ i or e s+ {Re t . 

Since m\N, this requires m\(e s — e s+ i + e t — £t+i)- Given the range of the 77/s and e/s, we have 
e s — e s+ i + e t — e t+ i = or ± m, which indicates that 

N N N N 
rj s - rj s+1 + 7] t - r] t+1 = 0, ±1, — — ± 1 , or ±1. (B.l) 

mm m m 

As discussed before, if we allow repeated elements in an m-pattern, the equivalence relation 1Z 
no longer necessitates a congruence relation as in pattern where each element is distinct. While the 
computation of high moments for general m-patterns appears intractable, fortunately we are able to 
illustrate how the difference in the modulo condition affects moment values by comparing the low 
moments for two simple patterns {a, b, a, b} and {a, a, 6, b}. 

B.2. The Fourth Moment. Although we will show that the higher moments differ by the way the 

elements are arranged in an m-pattern, the 4 th moment is in fact independent of the arrangement of 
elements. We first show that the 4 th moment for any m-pattern is determined solely by the frequency 
at which each element appears, and then show that this lemma fails for the 6 th moment and higher 
(see also lIXil ). Briefly, for the higher moments for patterns with repeated elements, there exist 
"obstructions to modulo equations" that make trivial some non-trivial configurations for patterns 
without repeated elements. Due to the obstructions to modulo equations, some configurations that 
are non-trivial for all-distinct patterns become trivial for patterns with repeated elements, making 
the higher moments for repeated patterns smaller. 

Lemma B.l. For an ensemble of real symmetric period m-block circulant matrices of size N, if 
within each m-pattern we have n i.i.d.r.v. {a r }" =1 , each of which has a fixed number of occurrences 
v r such that Y2r=i v r = m > 4 f/! moment of the limiting spectral distribution is 2 + Xlr=i( m) 3- 

By (12.101) . we calculate } +1 Y^i<i j k kn aijajkdkiaii for the 4 th moment. There are 2 ways of 

TV 2 

matching the 4 entries in 2 pairs: 

(1) (adjacent, 2 variations) a^- = a jk and a k i = an (or equivalently = an and a^ k = a k i); 

(2) (diagonal, 1 variation) = a k i and aj k = an. 

there are 3 matching s, with the two adjacent matching s contributing the same to the 4 th moment. 
We first consider one of the adjacent matchings, a^ = a jk and a k i = an. The pairing conditions 
(IB. 1 1) in this case are: 

(1) (diagonal condition) i — j = k — j (mod N), k — I = i — I (mod N); 

(2) (modulo condition) ilZk or jlZj, klZi or YR1. 
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Since 1 < i,j, k,l < N, the diagonal condition requires i = k, and then the modulo condition 
follows trivially, regardless of the m-pattern we study. Hence, we can choose j and / freely, each 
with N choices, i freely with N choices, and then k is fixed. This matching then contributes 
^ 3 = 1 (fully) to the 4 th moment, so does the other adjacent matching. 

We proceed to the diagonal matching, = clm and = an. The pairing conditions (IB. 1 1) in 
this case are: 

(1) (diagonal condition) i — j = I — k (mod N), j — k = i — I (mod N); 

(2) (modulo condition) iRX or jTZk, jTZi or kill. 

The diagonal condition j — k = i — I (mod N) is equivalent to % — j = I — k (mod N), which 
entails 

(1) i + k = j + 1, or 

(2) i + k = j + I + N, or 

(3) i + k = j + l-N. 

In any case, we only need to choose 3 indices out of i,j, I, k, and then the last one is fixed. In the 
following argument, without loss of generality, we choose I) and thus fix k. 

For a general m-pattern, we write % = Ar/i + e x , j = Arj 2 + e 2 , k = Ar/ 3 + e 3 , / = Arj 4 + e 4 , where 
V21 ^3) ^ {0, 1 , . . . , ^} and e\, e 2 , e 3 , e 4 E {0, 1, . . . , m — 1}. Before we consider the e/s, we 
note that there exist Diophantine constraints. For example, if i+k = j+l, given that 1 < i, j, I < N, 
k = j + I — i also needs to satisfy 1 < k < N. As a result, we need < r] 2 + Vi ~ Vi ^ x - 
Note that, due to the e/s, sometimes we may have < r] 2 + rj 4 — r]i < ^ + e, where the error 
term e G (— -j, y) and only trivially affects the number of choices of (772, 774, r/i) for a fixed m as 
N -+ 00. 

We now explore the Diophantine constraints for each variation of the diagonal condition (IB.2I) . 
The i + /c = j + 1 case is similar to that in HHMH . where, in a Toeplitz matrix, the diagonal condition 
only entails i + k = j + I, and there are obstructions to the system of Diophantine equations 
following the diagonal condition. However, the circulant structure that adds i + k= j + l + N and 
i + k = j + I — N to the diagonal condition fully makes up the Diophantine obstructions. This 
explains why the limiting spectral distribution for ensembles of circulant matrices has the moments 
of a Gaussian, while that for ensembles of Toeplitz matrices has smaller even moments. We now 
study the 3 possibilities of the diagonal condition for the circulant structure. 

(1) Consider i + k = j + I. We use Lemma 2.5 from HHMI to handle the obstructions to 
Diophantine equations, which says: Let I N = {1,...,N}. Then#{x,y,z e In ■ 1 < 
x + y-z<N} = |iV 3 + IN. 

In our case, let M = — . The number of possible combinations of (?7 2 , 774, 771) that allow 
< ?73 < ^ is |M 3 + |M0 For each of 772, 774, 771, we have m free choices of eg, and thus 
the number of is m 3 (§M 3 + \M) = lN 3 + 0(N). 

(2) Consider i + k = j + l + N. Note 1 < k < N requires 0<r/ 2 + ^4-^i + ^<^^-^< 
V2 + V-i ~ Vi — 0- Similar to the i + k = j + I case, we write M = — and S = 772 + r/4, and 
then < S - rji < ^ S < rn < M + S where obviously S < M. We have S + 1 
ways to choose (772, 774) s.t. 772 + 774 = S, and M — S + 1 choices of 77^ The number of 



In IIHMI . the related lemma is proven for 772 , ^74 , 771 <G N + , i.e., no cases where r\2 T \^ r \\ = 0- Thus we are supposed 
to start from 5 = 0; however, as N — > 00, the error from this becomes negligible. 
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(i, j, I) is thus 

M ( /If 3 ^ \ N 3 

m 3 ^(5 + l)(Af - 5 + 1) =m 3 I — + M 2 + ^M J = — + 0(N 2 ). (B.2) 

(3) Consider % + k = j + I — N. Now 1 < k < N requires < r/2 + 774 — r\\ — — < 
— =^ — < 772 + 774 — 771 < — . Again, we write M = — and S = r\\ + 774, and then 
M < S — r/i < 2M => S -2M < t]i < S - M where"obviously S > M. We have 
2M — S + 1 ways to choose (772, 774) s.t. 772 + 774 = S, and S — M + 1 choices of 771. The 
number of (i, j, I) is thus 

2M I ' M 3 \ N 3 

m 3 V (2M - S + IMS - M + 1) = m 3 + M 2 + -M ) = — + 0(N 2 ). (B.3) 

tti V 6 6 ) 6 

Therefore, with the additional diagonal conditions i + k = j + 1 + N and i + k = j + 1 — N induced 
by the circulant structure, the number of is of the order (| + | + l)^ 3 = N 3 , i.e. the 

circulant structure makes up the obstructions to Diophantine equations in the Toeplitz case. Since 
the ?7/s do not matter for the modulo condition, to make a non-trivial configuration, we may choose 
three 77/s freely, each with — choices, and then choose some q's that satisfy the modulo condition, 
which we will study below. 

For the modulo condition, it is necessary to figure out which zones the four entries are located in. 
Recall that the diagonal condition will always ensure that two paired entries are located in different 
areas. For the 4 th moment, each of the 3 variations of the diagonal condition is sufficient to ensure 
that any pair of entries involved are located in the right zones. We may check this rigorously by 
enumerating all possibilities of the zone-wise locations of the 4 entries, e.g. if i + k = j + / + N, 
if dij G zone 1, then a k i 6 zone 2^ As a result, for a pair of matrix elements in the diagonal 
matching, say = a ki , if i determines the slot in an m-pattern for and thus matters for the 
modulo condition, then I determines for a^; if j determines for a^, then k determines for a k i, and 
vice versa. 

With the zone-wise issues settled, we study how to obtain a non-trivial configuration for the 4 th 
moment. Recall the modulo condition for the diagonal matching: %RX or jlZk, jlZi or klZl. This 
entails 2 2 = 4 sets of equivalence relations, 

ininj;iTUTlk,jTlkni,jTlkni (B.4) 

Each set of equivalence relations appears with a certain probability, depending on the zone-wise 
locations of the 4 entries. For example, ilZllZj follows from ilZl and jTZi, which requires both 
and cijk E Area I. Regardless of the probability with which each set occurs, we choose one free 
index with N choices, and then another two indices such that these 3 indices are related to each 
other under 1Z. The number of choices of the two indices after the free one is determined solely by 
the number of occurrences of the elements in an m-pattern. 

We give a specific example of making a non-trivial configuration for the 4 th for two simple pat- 
terns {a, b, a, b} and {a, a, b, b}. Under the condition i + k = j + /, if 6 zone 1 and E zone 
3, then a^i E zone 4 and an E zone 2. We first select 771, r/2, 774 such that i, j, I and k = j + I — i 



This enumeration is complicated since the zone where an entry ay is located imposes restrictions on the choice of 

£ andj<f. 



e.g. when a, <G zone 2, we have i > ^ and j < ^ 
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k 


abab (theory) 


abab (observed) 


aabb (observed) 


abba (observed) 


N(0,1) 


2 


1.0000 


1.0016 


1.0014 


0.9972 


1 


4 


2.2500 


2.2583 


2.2541 


2.2405 


3 


6 


7.5000 


7.5577 


7.3212 


7.2938 


15 


8 


32.8125 


33.2506 


30.4822 


30.5631 


105 


10 


177.1880 


180.8270 


153.9530 


155.6930 


945 



Table 1 . Comparison of moments for various patterns involving a and b. The first 
column are the theoretical values for the moments of the pattern a, b, and the final are 
the moments of the standard normal. The middle three columns are 200 simulations 
of 4000 x 4000 matrices. 



satisfy the zone-wise locations^ In this case, based on pairing conditions (IB.ll) . pairing a^ = au 
and ajk = an will require ei7£e 4 and e 2 lZei, or equivalently ei7^e 2 7?.e4. Without loss of generality, 
we can start with a free e\ with 4 choices, then there are 2 free choices for each of e 2 and e 4 , and 
then we have a non-trivial configuration. We have similar stories under the other two variations of 
the diagonal condition and with other zone- wise locations of and a ki . Therefore, we can choose 
three out of four 77/ s freely, each with ^ choices, then one q with 4 choices, then another two 
e/s each with 2 choices, and finally the last index is determined under the diagonal condition. As 
discussed before, such a choice of indices will always satisfy the zone-wise requirements and thus 
the e-based pairing conditions. Thus there are (^) 3 ■ 4 • 2 • 2 = ^- choices of k, I) that will 
produce a non-trivial configuration. It follows that the contribution from the diagonal matching to 
the4*rnomentis^(f + ± + |)f = \. 

The computation of the 4 th moment for the simple patterns {a, b, a, b} and {a, a, b, b} can be 
immediately generalized to the 4 th moment for other patterns. As emphasized before, both adjacent 
matchings contribute fully to the 4 th moment regardless of the m-pattern. For diagonal matching, 
the system of Diophantine equations induced by the diagonal condition are also independent of the 
m-pattern in question, and the way we count possible configurations can be easily generalized to an 
arbitary m-pattern. We have thus proved Lemma |BTl 

Note that Lemma IB.ll implies that the 4 th moment for any pattern depends solely on the fre- 
quency at which each element appears in an m-period. Besides the {a, a, b, b} pattern that we have 
studied in depth, we may easily test two extreme cases. One case where n = m, i.e. each ran- 
dom variable appears only once, represents the m-block circulant matrices from Theorem 11.41 for 
which the 4 th moment is 2 + (and m — 1 represents the circulant matrices for which the 4 th mo- 
ment is 3). Numerical simulations for numerous patterns including {a, a, b}, {a, b, b}, {a, b, b, a], 
{a, b, c, a, b, c}, {a, b, c, d, e, e, d, c, b, a} et cetera support Lemma |BTI as well; we present results of 
some simulations in Tables Q] to [3] 



It is noteworthy that the specific location of an element still depends on the e/s, but as N — > oo, the probability 
that the 77/s alone determine the zone-wise locations of elements approaches 1, i.e. the probability that adding the e/s 
changes the zone-wise location of an element approaches 0. 
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k 


ababab 


aaabbb 


aaaabbbb 


aaaaabbbbb 


aababb 


2 


1.0000 


1.0008 


1.0001 


0.9984 


0.9996 


4 


2.2500 


2.2541 


2.2441 


2.2449 


2.2502 


6 


7.5000 


7.3011 


7.2098 


7.2551 


7.2319 


8 


32.8125 


30.3744 


29.5004 


30.0127 


29.5378 


10 


177.1880 


155.0380 


145.8240 


150.7220 


145.4910 



Table 2. Comparison of moments for various patterns involving a and b. The first 
column are the theoretical values for the moments of the pattern a, b. The remaining 
columns are 200 simulations of 3600 x 3600 matrices. 



k 


abcabc 


abccba 


aabbcc 


abbcca 


aabcbc 


2 


1.0000 


1.0005 


1.0006 


0.9983 


1.0013 


4 


2.1111 


2.1122 


2.1153 


2.1047 


2.1161 


6 


6.1111 


6.0248 


6.0540 


6.0083 


6.0235 


8 


22.0370 


20.9398 


21.2004 


20.9908 


20.8411 


10 


94.6296 


85.0241 


87.0857 


85.9902 


84.2097 



Table 3. Comparison of moments for various patterns involving a and b. The 
first column are the theoretical values for the moments of the pattern a, b, c. The 
remaining columns are 200 simulations of 3600 x 3600 matrices. 



B.3. The Sixth Moment. Although for an m-pattern with each element appearing at a fixed fre- 
quency, the 4 th moment is independent of how the elements are arranged within the pattern, the way 
the elements are arranged in an m-pattern does affect higher moments and thus the limiting spectral 
distribution. As we will show for the 6 th moment, for patterns with repeated elements, there exist 
"obstructions to modulo equations" that make trivial some non-trivial configurations for patterns 
without repeated elements. We illustrate this by explicitly computing the 6 th moment for the pattern 
{a, b, a, b}, and then showing why the 6 th moment for {a, a, b, b} differs. It will then be clear that 
the modulo obstructions persist for more complicated patterns and higher moments. 

For the 6 th moment, we calculate J2i<i,j,k,i,m,n<N a ij a jkakiai m amna n i by (12.101) . There are 
(6 — 1)!! = 15 matchings, which can be classified into 5 types, so that the 6 entries are matched in 
3 pairs: 

(1) a i:j = a jk , a u = a lm , a mn = a ni (adjacent, 2 variations). 

(2) a i:j = a jk , a u = a ni , a lm = a mn (semi-adjacent- 1, 3 variations). 

(3) ay = a jk , a ki = a mn , a tm = a ni (semi-adjaent-2, 6 variations). 

(4) ciij = ai m , a jk = a ni , a kl = a mn (diagonal- 1, 3 variations). 

(5) ciij = ai m , a jk = a mn , a ki = a ni (diagonal-2, 1 variation). 

Similar to the 4 th moment computation, we first take advantage of adjacent cases. For example, 
if = aj k , then the two pairing conditions ( IB.ll) require 

(1) i — j = k — j i = k. Given that i, j, k, I 6 {1, 2, . . . , N}, neither %— j — k — j + N nor 
i — j = k — j — N is possible. 

(2) % = k (mod 2) or j = j (mod 2), depending on the zone- wise location of and dj k . 
Either follows trivially from the previous condition. 
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For Type 1 (adjacent), take = aj k , a k i = ai m , a mn = a ni , the diagonal condition requires 

i — j — k — j,k — I — m — l,m — n — i — n =>- i — m — k. (B.5) 

By the discussion on the adjacent case, the modulo condition is satisfied trivially. We can then 
freely choose i, j, I, n, each with N choices, and make a non-trivial configuration that contributes 
^ 4 = 1 (fully). Type 1 matchings thus contribute 2 x 1 = 2 (2 variations of Type 1) to the 6 th 
moment. 

For Type 2 (semi-adjacent- 1), take a^ = aj k , a k i = a n i, ai m = a mn , the adjacent case = aj k 
requires i = k as discussed before. Thus the second pair a k i = a ni is equivalent to a k i = a nk , which 
is again an adjacent case. The third pair ai m = a mn is an adjacent case itself. Thus Type 2 is in fact 
equivalent to Type 1, and contributes 3 x 1 = 3 to the 6 th moment. 

For Type 3 (semi-adjacent-2), the adjacent case = aj k requires i = k as discussed before. 
Thus the third pair ai m = a ni is equivalent to ai m = a nk , and the second and the third pair combined 
make the diagonal matching as in the 4 th moment computation. We have shown that this diagonal 
matching contributes | to the 4 th moment (see Lemma IBTTT) . Note that j is free with N choices 
despite the restriction i = k. Thus this matching contributes |, and this type 6 X | = |, to the 6 th 
moment. 

Note that Type 1 and 2 are independent of the m-block circulant pattern along the diagonals 
in an m-block circulant matrix, and Type 3 also applies to other variations of {a, b, a, b} such as 
{a, a, b, b} and {a, b, b, a}. Type 1 through 3 combined, we have 2 + 3 + | = 6.5 in the 6 th moment. 

We proceed the diagonal matchings, for which we will discuss the modulo obstructions, and start 
with a simple case for Type 4. Take the matching = ai m , ajk = a ni , an = a mn as an example, 
the two pairing conditions (IB.ll) require: 

(1) i — j = m — I, j — k — i — n, k — I = n — m i — j — m — I = n — &0 which shows 
that we need to choose only 4 of the 6 indices, and the other 2 are determined. 

(2) ilZm or jTZl, jlZi or klZn, klZn or llZm, depending on the zone-wise locations of the 
entries. For example, if E zone 1, then % — j = m — Z =>- a/ m e zone 4. We have 2 3 = 8 
sets of equivalance relations, categorized as follows. Category (1) (4 sets): ilZmlZj, klZn; 
jlZllZm, klZn; ilZmlZl, klZn; jlZllZi, klZn. 

Category (2) (2 sets): ilZmlZjTU; jTUIZilZm. 
Category (3) (2 sets): ilZm, klZn; jTZl, klZn. 

Each set of equivalance relations appears with a certain probability, and the probabilities of ob- 
serving each 1Z set sum up to 1. We show below that, regardless of the probability of observing 
each set, each set contributes | to the 6 th moment, and thus the probability-weighted contribution is 
simply \. 

For Cat.(l), the set of equivalence relations ilZmlZj, klZn requires a^-, ajk, a k i G zone 1 or 3. 
Thus we can start with a free i with N choices, then select m,j, each with y choices, such that 
ilZjlZm. Then we pick a k, and note that i — j = n — k, ilZj =^ klZn. Recall that, for the pattern 
{a, b, a, b}, ilZj indicates 2\{i — j), and it follows that 2\{n — k). In other words, we can freely 
choose a k with N choices, and the diagonal condition i — j = n — k ensures that we have a good 



We temporarily ignore i — j = m — I + N and i — j = m — I — N for simplicity. In fact, as we show in the 4 
moment computation, the i — j ~ m — I + N case and the i — j = m — I + N case, each of which has ^- + 0(N 2 ) 
solutions, together make up the obstructions to a Diophantine equation like i — j = m — I that has + 0{N 2 ) 
solutions. 



34 



MURAT KOLOGLU, GENE S. KOPP, AND STEVEN J. MILLER 



n. This set thus contributes 4j • (N • y ■ y • N) = \. The same analysis applies to the other 3 sets 
in Cat.(l). 

For Cat.(2), take the set iR.rnR.jJU . We start with a free i, and then select m,j, each with — 
choices, such that iRrriRj. Note that, again, % — j = m — /, iRj =>■ mJU iRmlZjJU. This set 
thus contributes -jb • (N ■ y • y • TV) = \. The same analysis applies to the other set in Cat. (2). 

For Cat. (3), take the set ilZm, klZn. We start with a free i, and then select m with — free choices 
such that iRm. Then we choose a free k with N choices and n with y choices such that klZn. This 
set thus contributes • (N ■ y • N • y) = \. The same analysis applies to the other set of Cat. (3). 

Since each individual set of equivalence relations in Cat.(l)-(3) contributes equally, the probability- 
weighted contribution to the 6 th moment is |. Therefore, Type 4, with 3 variations, contributes |. 

Similarly, the pairing conditions (IB.ll) entail the following for Type 5 (diagoal 2). 

(1) i — j = m — I = k — n. 

(2) 2 categories of equivalence relation set. 

Cat.(l)(6 sets): ilZmlZk, jlZn; jJUlZn, klZm; ilZm, jlZnJU; jJU, kJZmJU; iJLmJtk, llZn; 
jJUlZn, klZi. 

Cat.(2)(2 sets): iJlmJlk; jJUJln. 

Replicating the analysis of Type 4, we find that, since each set of equivalence relations in Cat.(l) 
and Cat. (2) contributes ~, the probability-weighted contribution must be | as well. Since Type 5 
has only 1 variation, Type 5 contributes \ to the 6 th moment. 

Therefore, the combined contribution from Type 4 and Type 5 is | + \ — 1. The 6 th moment for 
the pattern {a, b, a, b} is then 6.5 + 1 = 7.5. 

Now we examine why the contribution from diagonal matchings for the pattern {a, a, b, b} differs 
from that for {a, b, a, b}. As discussed before, Type 1 through 3 matchings, with a total contribution 
of 6.5, also apply to {a, a, b, b}. For Type 4 and 5, however, the combined contribution is less than 
1. Recall a key argument in the analysis of Type 4 matching before: for the 2-block circulant 
{a, b, a, b} pattern, under i — j = m — I = n — k, if we choose ilZjlZm, i.e. i = j = m (mod 2), 
then / = j+m—i will satisfy / = m (mod 2) as well. Namely, %—j = m—l, %Rj =^ mJU. However, 
for {a, a, b, b}, if we specify 1Z as sJlt <^=^ mod (s,4), mod (t,4) e {1,2} or mod (s,4), 
mod (t, 4) G {0, 3}, and choose ilZjlZm, it is possible that / = j + m — i is not related to m 
under 1Z. For instance, when mod (i,4) = 2, mod (j, 4) = 1, mod (m,4) = 3, we have 
iJlj, but mod (Z,4) = 2. Some configurations that are non-trivial for {a, b, a, b} then become 
trivial for {a, a, b, b}, while all the non-trivial configurations for {a, a, b, b} are still non-trivial for 
{a, b, a, b}. Thus, we expect the 6 th moment for {a, a, b, b} to be smaller than that for {a, b, a, b}, 
which is also evidenced by numerics. We phrase such a loss of non-trivial configurations as due to 
"obstructions to modulo equations", or "modulo obstructions" for short, which will clearly persist 
in higher moments for general m-block circulant patterns with repeated elements. 

Based on the brute-force computation above, we may also bound the even moments for gener- 
alized m-block circulant patterns. It is clear that a lower bound is the moment for the m-block 
circulant pattern of the same period length and in which each element is distinct. For example, in 
terms of high (2k th , k > 2) moments, {a, a, b, b} > {a, b, c, d} (both of length 4). In the compu- 
tation of high moments, a pattern with repeated elements has all the non-trivial configurations that 
an all-distinct pattern of the same length can have, and gains extra non-trivial configurations due to 
the repeated elements. 

An easy upper bound is the moment of the standard Gaussian, which is the limiting spectral dis- 
tribution for the ensemble of circulant matrices. We may also easily find a sharper upper bound 
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for a family of simple m-block circulant patterns in which each element appears at the same fre- 
quency, e.g. {a, b, c, c, b, a}, {a, a, b, c, b, c}, etc. For this family, an upper bound will be associated 
with a pattern where each element only appears once. For example, in terms of high moments, 
{a, b, c} > {a, b, c, c, 6, a}. We may take {a, b, c} as {a, 6, c, a, b, c}, and note that, although in 
{a, 6, c, a, 6, c}, the probability of choosing each letter is the same as in {a, b, c, c, b, a}, the former 
pattern is free of modulo obstructions that exist for the latter. 

For a more general m-block circulant pattern, however, a sharper upper bound is not easily attain- 
able. For instance, it is not clear whether {a, a, b, c] > {a, b, c}. Some numeric evidence suggests 
that a pattern in which gcd(vi, v 2 . . . vg) = 1, where vg is the number of occurrences of an element 
in an m-period, has larger high moments than those with the same frequency of each element but 
gcd(ui, v 2 . . . vi) > 2. For example, {a, b, c, c} > {a, a, b, b, c, c, c, c}. 

Obviously, the accounting above will become significantly more involved for more complicated 
patterns or higher moments, but the basic ideas remain the same. We also foresee that as the 
moments get higher, the number of configurations that contribute trivially will increase so quickly 
that the higher moments get increasingly farther below the standard Gaussian moments. This is also 
evidenced by simulations. 

B.4. Existence and Convergence of High Moments. Although it is impractical to find every mo- 
ment for a general m-block circulant pattern using brute-force computation, we are still able to 
prove that, for any m-block circulant pattern, every moment exists, is finite (and satisfies certain 
bounds), and that there exists a limiting spectral distribution. In addition, the empirical spectral 
measure of a typical real symmetric m-block circulant matrix converge to this limiting measure, 
and we have convergence in probability and almost sure convergence. 

We have shown that all the odd moments vanish as N — > oo, and thus we focus on the even 
moments. We need to prove the following theorem. 

Theorem B.2. For any patterned m-block circulant matrix ensemble, liniTv^oo M 2 k(N) exists and 
is finite. 

Proof. It is trivial that M 2 k(N) is finite. As discussed before, it is bounded below by the 2fc th 
moment for the ensemble of m-block circulant matrices where, in the m-pattern, each element is 
distinct, and more importantly it is bounded above by the 2k th moment for the ensemble of circulant 
matrices, and we know that the limiting spectral distribution for this matrix ensemble is a Gaussian. 

We now show that liniAr^oo M 2 k(N) exists. To calculate M 2 k(N), we match 2k elements from 
the matrix, {a^, di 2 i 3 , . . . , a^,,^}, in k pairs, where it e {1, 2, . . . , N} and this will give (2k — 1)!! 
matchings. For each matching, there are a certain number of configurations, and most of such 
configurations do not contribute to the moments as N — > oo. 

For the m-block circulant pattern, the equivalence relation 1Z implies that e s lZe t+ i e s = e t+ i, 
and since m| (e s — e s+1 + e t — e t+ i), we have e s+1 = e t as well (see (1B.1| ))F1 Thus r] s — i] s+ i + rj t — 
7] t +i = or ± — , three equations that have (^) 3 + 0((^) 2 ) solutions in total, as we have shown 
in the 4 th moment computation. 

However, if there are repeated elements in an m-period, then e s 7Ze t +i no longer necessitates 
e s = e t+1 , and it is possible that (e s — e s+ i + e t — e t+1 ) = ±m. Thus, the zone-wise locations of 
elements matter in making non-trivial configurations. Recall that the zone-wise location (see (IB.ll) ) 
of an element determined by (i s +i — i s )' if a i s i s+1 is m zone 1 or 3 (Area I), e s determines 

1 1 This explains why, for an m-pattern without repeated elements, the zone-wise locations of matrix entries do not 
matter in making a non-trivial configuration. 
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the slot of a isis+1 in an m-period; if a isis+1 is in zone 2 or 4 (Area II), e s+1 determines the slot of 
a isis+1 in an m-period. In addition, the diagonal condition will always ensure that two paired entries 
a isia+1 and a itit+1 are located in different areas. 

Recall that for any matching M. , the k pairs of matrix elements, each pair in the form of cbi 3 i s+1 = 
a>nh+i> are fixed. For any Ai, to make a non-trivial configuration, we first choose an e vector of 
length 2k. If we choose all the e/s freely, there are m 2k possible choices for an e vector, most of 
which do not meet the modulo condition, and trivially, m 2k is an upper bound for the number of 
valid e vectors. It is noteworthy that out of the 2k q's of an e vector, only some of the q's will matter 
for the modulo condition. Which q's in fact matter depends on how we pair the 2k matrix entries 
Oisis+i's an d the zone-wise locations of the paired aj s j s+1 's, which we cannot determine without 
fixing the 77/s (and thus the i/s). 

However, for any matching, the way we pair the 2k matrix entries into k pairs is fixed, and for 
each fixed pair a isis+1 = a itit+1 , two q's will matter for the modulo condition: either e s 1Ze t+ i or 
e s+ {R,e t . Thus there are 2 k ways to choose k pairs of q's for each matching. For each way of 
fixing the k pairs of e/s, we examine each e pair, say (e^, q 2 ), and there are a certain number of 
choices of (e^, q 2 ) such that e^TZe^. Continuing in this way, for each e pair, we choose two e/s 
that satisfy the equivalence relation 1Z. Note that an ei may matter twice, once, or never for the 
modulo condition depending on the zone- wise locations of the Oi a i i+1 's. We then choose the other 
e/s that do not matter for the modulo condition such that for each pair of di s i s+1 = Oi t i t+1 , we have 
e s — e s+ i + e t — e t +\ = or ± m, and finally we have a valid e vector. The number of valid e vectors 
will be determined by m, k, and the pattern of an m-period, but will be independent of N since the 
system of k equivalence relations for the modulo condition does not involve N. 

With a valid e vector, we have fixed the zone-wise locations of the 2k matrix elements by fixing 
the q's that matter for the modulo condition. We now turn to the diagonal condition and study the 
rjiS. With k equations in the form of 

m(r] s - 7] s+ i +r) t - rjt+x) + (e s - e s+ i + e t - e t+1 ) = or ±N, (B.6) 

and (e s — e s+1 + e t — e t+ i) known in each of the k equations, we in fact have k equations in the form 
of 

7] s - 7] s+ i +7] t - Vt+i = 7, (B.7) 

where 7 £ {0, ±1, — , — ± 1, — — , — — ± 1}. This gives us k + 1 degrees of freedom in choosing 
the 77/ s, and trivially, we can have at most vectors of 77/s. Since the e vector is fixed, for 

one equation 7] s — 7? s+1 + r\ t — rj t +i = 7, there are only 3 choices of 7. With k equations in this 
form, we have at most 3 fc systems of 77 equations. Note that not all of the 77 vectors satisfying an 
7] equation system derived from the diagonal condition will help make a non-trivial configuration, 
since the 77/s need to be chosen such that the resulted aj s j s+1 's will satisfy the zone- wise locations 
in order to be coherent with the pre-determined e vector. For example, if in a pair of matrix entries 
a-isis+i = a ith+i where e s TZe t+ i, even though the 77/s are chosen such that 7? s - r] s+ i + r] t - rj t+1 = 7, 
it is possible that a isis+1 , a itit+1 are located in certain zones such that we need e s+ {Re t to ensure a 
non-trivial configuration. 

The following steps mirror those in HHML Denote an 77 equation system by S. For any S we 
have k equations with 771, 772, ... , r] 2k £ {1, 2, . . . , Let z t = ^ £ ^, . . . , 1}. Without 
the zone-wise concerns discussed before, the system of k equations would have k + 1 degrees of 
freedom and determine a nice region in the (k + 1) -dimensional unit cube. Taking into account 
the zone-wise concerns, however, we will still have k + 1 degrees of freedom. For example, for 
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a pair of matrix elements a isiB+1 = a itit+1 , the system S requires 77*, — r] s+ i + rj t — rjt+i = 7. If 
we need e s lZe t+ i to make a non-trivial configuration, say a isia+l E zone 1, then we will obtain 
an additional equation < i s+ i — i s < y — 1 =^ < (r) s+ i — r) s ) + e s+i — e s < y — 1 with 
(e s+ i — e s ) G { — m + 1, — m + 2, . . . , 0, 1, . . . , m — 2, m — 1}. Based on the region determined 
by rj s — 7] s+ i + 7) t — r] t+1 = 7, this additional zone-related restriction will only allow a slice of the 
region for us to choose valid 77/s. With k zone-wise restrictions, only a proportion of the original 
region in the unit cube will be preserved for the choice of the 77 vector. Nevertheless, the "width" of 
each slice is of order y, and we still have k + 1 degrees of freedom. 

Therefore, with m fixed and as 7Y — > 00, we obtain to first order the volume of this region, 
which is finite. Unfolding back to the r//s, we obtain M 2k (S)(^) k+1 + O fc ((^) fc ), where M 2k (S) 
is the volume associated with this 77 system. Summing over all rj systems, we obtain the number 
of non-trivial configurations for the 2k th moment from this particular e vector. Next, within a given 
matching Ai, we sum over all valid e vectors, the number of which is independent of iV as we have 
shown before. In the end, we sum over the (2k - 1)!! matchings to obtain M 2k N k+1 + O k (N k ), and 
the 2k th moment is simply M2fcjV ^I+° fc(Affc) = M 2k + O(^). □ 

The above proves the existence of the moments. The convergence proof follows with only minor 
changes to the convergence proofs from [|HM[|MMSl . 
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